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Abstract 



A neutrinophilic Higgs model has tiny vacuum expectation value (VEV), 
which can naturally explain tiny masses of neutrinos. There is a large 
^ energy scale hierarchy between a VEV of the neutrinophilic Higgs doublet 

and that of usual standard model-like Higgs doublet. In this paper we at 
^ first analyze vacuum structures of Higgs potential in both supersymmetry 

(-r^ (SUSY) and non-SUSY neutrinophilic Higgs models, and next investigate 

O a stability of this VEV hierarchy against radiative corrections. We will 

S^ show that the VEV hierarchy is stable against radiative corrections in 

[^ both Dirac neutrino and Majorana neutrino scenarios in both SUSY and 

O non-SUSY neutrinophilic Higgs doublet models. 



X 



1 Introduction 

The recent neutrino oscillation experiments gradually reveal a structure of lepton sector [T],|2]. 
However, from the theoretical point of view, smallness of neutrino mass is still a mystery and 
it is one of the most important clues to find new physics beyond the standard model (SM). 
A lot of ideas have been suggested to explain the smallness of neutrino masses comparing to 
those of quarks and charged leptons. How about considering a possibility that the smallness 
of the neutrino masses is originating from an extra Higgs doublet with a tiny vacuum expec- 
tation value (VEV). This idea is that neutrino masses are much smaller than other fermions 
because the origin of them comes from different VEV of different Higgs doublet, and then 
we do not need extremely tiny neutrino Yukawa coupling constants. This kind of model is 
so-called neutrinophilic Higgs doublet model [3]-[l3], where a neutrinophilic Higgs take a VEV 
of (9(0.1) eV in Dirac neutrino scenario[5l El El |9], while a VEV of 0{1) MeV in Majorana 
neutrino scenario with TeV-scale seesaw[3l IH [71 UHl El [El [13]. The non-supersymmetric 
(non-SUSY) neutrinophilic Higgs doublet model is sometimes called z/THDM. The (collider) 
phenomenology in z/THDM is interesting, since a charged Higgs boson is almost originated 
from the extra neutrinophilic Higgs doublet and its couplings to neutrinos are not small. The 
characteristic signals of the z/THDM could be detected at LHC and ILC experiments [HI [TT] . 
Not small neutrino Yukawa couplings in the z/THDM can also make low energy thermal lepto- 
genesis work^l2j. The SUSY version of neutrinophilic Higgs doublet model have been suggested 
in Refs. fTH |T3], where a thermal leptogenesis in a low energy scale works without gravitino 
problem [la [JSQ 

Anyhow, a neutrinophilic Higgs model has tiny VEV, and there is a large energy scale 
hierarchy between a VEV of the neutrinophilic Higgs doublet and that of usual SM-like Higgs 
doublet. In this paper, we at first analyze vacuum structures of Higgs potential in both 
SUSY and non-SUSY neutrinophilic Higgs models, and next investigate a stability of this 
VEV hierarchy against radiative corrections. We will show that the VEV hierarchy is stable 
against radiative corrections in both Dirac neutrino and Majorana neutrino scenarios in both 
SUSY and non-SUSY neutrinophilic Higgs doublet models. 

2 z/THDM 

Let us analyze vacuum structures of Higgs potential in non-SUSY neutrinophilic Higgs model, 
i.e., z/THDM at first, and next investigate a stability of this VEV hierarchy against radiative 
corrections. 

* Cosmological constraints were argued in Ref.jl^, however, a setup of them is different from usual neu- 
trinophilic Higgs doublet models, since it includes a light Higgs particle. 



2.1 Vacuum structure in tree-level potential 



We here overview the i/THDM, where we introduce a neutrinophihc Higgs doublet $;^ and 
Z2-parity as follows. 



Fields 



SM Higgs $ 
neutrinophihc Higgs $1, 
right-handed neutrino 
others 



Z2-parity 



+ 



Lepton number 






1 
± 1: leptons, 0: quarks 



Yukawa interactions are given by 

-Cyukawa = V'^Ql^Ur + V^Ql^Dr + v'L^Er + y^L^.N + h.C. 



(2.1) 



where $ = 2(72$, and generation indexes are omitted. Note that the right-handed neutrino 
only couples with $1, through the Yukawa interaction, and this is the origin of smallness of 
the neutrino masses. When we include Majorana mass of right-handed neutrinos ^MN'^N, 



this model becomes Majorana neutrino scenario through the seesaw mechanism 
potential is given by 



A Higgs 



ri^THDM 



yurnu., = -m2$t$ + m^$t$, - m2($t$, + $t $) + ^($t$)2 + ^($t$,)2 

A5 



+ A3($t<l')($|:<l'.) + A4($t<l'.)($t'^) + y [('^^'^.)' + ('^t'^)'], (2.2) 

where parameters are mi ~ m2 ~ O(IOO) GeV and Aj ~ C>{1) (i = 1,...,5). As for a 
magnitude of 1^.3 1, we take {O{10'^'^) GeV)^ for Majorana neutrino scenario, and {O{10~^) 
MeV)^ for Dirac neutrino scenario. Notice that $ has negative mass squared as {—mf) < 0. 
The Higgs doublets are assumed to be take real VEVs as ($) = (t'i,0)"^ and ($j,) = (f2,0)"^, 
then, stationary conditions are given by 











1 (}\/i'THDM 

2 dvi ~ 
1 dV'^'^^^^ 

2 d^2 



-mlvi — m\v2 + Aif^ + Afifg, 
-- ■m\v2 - m\vi + A2t;2 + ^'vlv2-, 



(2.3) 
(2.4) 



where A = A3 + A4 + A5. We sort the following three cases by magnitude relations between 
|fi| and |f2|. 

1- \v2\ ^ \vi\ case: This vacuum is what the i/THDM wants to realize. The magnitudes of 
VEVs are given by 

hl^W-1, t;2^ /^ , (2.5) 

V '^1 m2 + Xvf 



and a potential height at the vacuum is given by 

T/THDM 2 2 , ^1 4 ^ ^' 



2Ai 



2. |fi| <ti \v2\ case: This vacuum suggests V2{m2 + A2f|) = from Eq.(2.3), and thus 



2 ^ ^ K > 0) 

X? K<o) 



^2 = S m2 , o _, (2.6) 



^o / 2 



The case of f ^ = contradicts \vi\ <^ |f2|. Another case of fg = — m2/A2 suggests the 
value of f 1 as vf = m|f |/(— mf + Af|)^ > 0, where a potential height is given by 



Ai^THDM "^'' 



^\v2\Mvi\ - 2A ■ ^ '^ 



3. I^il ~ |f2| case: Neglecting tiny parameter m^, the stationary conditions Eqs.(2.3) and 



(2.4) become 



mlvi + Xiv'l + Xvivl = 0, (2.8) 

mlv2 + Asi^a + ^^^i^2 = 0. (2.9) 



Then, VEVs are given by 



^ _A2m[+Am| ^ Am[+Aim| 

A2-A1A2 A2-A1A2 



and the potential height at the vacuum is estimated as 



ruTUBM „. Aim^ + A2m^ + 2\m \ml 
2(A2 - A1A2) 



yvTimU ^ /\l//62 -r /V2"n T ^/v»ti»62 (2 111 



\V2\ < |t;i| should be a global minimum. Conditions of V^J^^^^^ < V;:;™Jf or V^^^^^^^ < 
*^bih|i'2| suggesx 

P = (A3 + A4 + A5)'>AiA2. (2.12) 

This is a necessary condition for vi ^ f 2 to be the global minimum, and an additional condition 

4 4 

— ^ < — ^ makes the vacuum true global minimum. For the potential to be bounded from 
below [TOl [16], quartic terms must satisfy 

x/A^>-(A3 + A4±A5), V'aIA^>-A3, Ai>0, A2>0. (2.13) 



These are the conditions of bounded below of the Higgs potential. We can show that a case 
of A < cannot satisfy the global minimum condition. Therefore, only a case of A > can 
satisfy the global minimum condition. Thus, in order for the desirable vacuum vi ^ f 2 to be 
the global minimum, a condition 

0<v^A^<A, V'aIa^>-(A3 + A4-A5), Ai,A2>0 (2.14) 

is needed. 

Next, let us estimate a curvature (mass squared) at each vacuum, which is given by 

^2 _ 1 gy™"^ _ f-nil + 3Ait;2 + \vl -ml + 2\v^V2 
'^ ~ 2 dvidvj ~ V -ml + 2At;it;2 m| + 3A2t;i + Xvf 



Then, the eigenvalue equation (eigenvalue: x) is given by 

x^ - {-ml + ml + (3Ai + X)vl + (3A2 + X)vl)x - mjml - m^ + 3A(Ait;^ + A2t;2) 

+ (3Aim^ - Xml)vl - (3A2m? - Am^)f^ + 3(3AiA2 - X^)vlvi + 4m2AAiA2 = 0, (2.16) 

and we can estimate the curvature for above three cases. 



1. 1^1 1 3> \v2\ case: The eigenvalue equation in Eq.(2.16) becomes 

x^ - {-ml + ml + (3Ai + X)vl)x + 3AAif^ + (3Aim^ - Xml)vl - mlm\ ~ 0. (2.17) 



By using Eq.(2.5), it becomes 

{x - 2ml){x - Xvl + ml) = 0, (2. 18) 

which means 

x = 2ml, Xvl + mj. (2.19) 

Thus, m^l > and Amf + Aim^ must be needed for x > 0. 



2. |t>i| <^ |t>2| case: Using V2 = —^ in Eq.(2.6), the eigenvalue equation in Eq.(2.16) 
becomes 

{x + 2ml){x-{Xvl-ml)) = 0. (2.20) 

Thus, the solution is given by 

x = -2ml, Xvl -ml, (2.21) 

which means m| < 0, Am^ + A2mf < for x > 0. 



3. \vi\ ~ |t>2| case: By neglecting mo and using Eq.(2.10), the eigenvalue equation in 



Eq.(2.16) becomes 

x^ - 2(Aii;2 + \2vl)x - 4(P - \i\2)vlvl = 0, (2.22) 

which means two eigenvalues xi,X2 should satisfy 

xi + X2 = 2(Ait;i2 + A2W2), xiX2 = -4:{\^ - \i\2)vlvl (2.23) 

Since positive xi,X2 give positive xi + X2,xiX2, the positive curvature condition at this 
vacuum is given by 

Xivl + X2VI > 0, (2.24) 

-(P - \i\2)vfvl > 0. (2.25) 

Thus, A^ — A1A2 < 0, is a positive curvature condition at the vacuum of |fi| ~ \v2\. 

The squared masses of the charged Higgs and of the pseudoscalar must be also positive. These 
conditions are equivalent to 

ml + \2vl + X^vl > (2.26) 

ml + \2vl + (A3 + A4 - X^Wi > 0. (2.27) 

Summarizing conditions for the vacuum we want, at first, A^ — A2A2 must be positive for 

4 4 

the vacua of |fi| ^ \v2\ and |fi| ^ \v2\ to be lower than that of |fi| ~ \v2\, and ~^ < ~^ 
makes the vacuum of |i;i| ^ \v2\ the global minimum. Note that A must be also positive to be 
consistent with the conditions of the potential bounded from below. Next, positive curvature 
conditions are mil > or Xml + Aim2 > with mil < 0. Finally, positive curvature of the 
charged Higgs and the pseudoscalar components require m^ + Xsvf > and m^ + (A3 + A4 — 
As)^^ > at |fi| ^ \v2\- In Table 1, we show which vacuum becomes the global minimum 
depending on signs of m^. A, and A^ — A1A2. 

Can a "local minimum" at \v2\ <^ |fi| in (2), (4) and (5) be our vacuum? It might be 
possible if a life time of the local minimum is long enough. There is a transition process from 
the local minimum at \v2\ -C |fi| to the global minimum at |fi| ~ \v2\- Its transition probability 
of tunneling rate suggests the life time is much shorter than an age of our universe, since a 
"distance" and a "height" of wall between the local and global minimums are both (9(100) 
GeV with 0{1) couplings of Aj in Higgs potential. So, unfortunately, the local minimum 

cannot be our vacuum. Therefore, in the i/THDM, we must use the suitable parameter setup 

4 4 

as (1) or (3) with -^ < -^. 

Before closing this subsection, we comment on recent analyzes of vacuum structure in 

general THDM. For example, in Ref.[T6], they investigated the vacuum instability of charge 

and/or CP breakings at tree level. As for so-called Inert Doublet Model (IDM) [17], it has exact 

Z2-symmetry with r/ig = 0. This Inert Doublet does not couple with any matter fermions, 

which is crucial difference from our model. 
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Table 1: Six cases which satisfy conditions in Eq.(2.13), (2.26) and (2.27). Here GM means 



"Global Minimum" and PC means "Positive Curvature", and / in GM (PC) says each vacuum 
can be the global minimum (has positive curvature), (a) and (b) mean that the positive 
curvature requires conditions of (a): Xml + \iml > and (b): —Xm"^ — \2'm,l > 0, respectively. 
Two simultaneous / in GM means |fi| ^ \v2\ (|fi| ^ |f2|) vacuum becomes the global 

A A A A 

minimum when 



2Ai 



< 



2A2 



2Ai 



> 



2A2' 



2.2 Stability against radiative corrections 

Now we are in a position to investigate the stability of the VEV hierarchy \v2\ ^ \vi\ against 
radiative corrections. First of all, we should remind that the small magnitude of \'mW plays a 
crucial role for generating the tiny VEV of |f2|(^ l"^!!)- Its smallness is guaranteed against 
radiative corrections, since it is the "soft" breaking mass parameter of the Z2-symmetry. As 
noted in Ref.[8], the radiative correction to this parameter is expected to be logarithmic. 
For analyses of the vacuum stability, we should use Coleman- Weinberg type 1-loop effective 
potential[T8], and analyze the stability of the VEV hierarchy. This 1-loop effective potential 
contains infinite number of irrelevant operators with zero-momentum Higgs fields in the exter- 
nal lines, and is calculated by a summation of them. However, for the investigation of stability 
of the VEV hierarchy, it is enough for us to pick up only diagrams which have external lines 
of mixture of $ and $j,. Furthermore, we should notice that, when one ^^ is added in the 
external lines, a coefficient of the effective operator should have suppression factor, |f2/'n^i,2|- 
Thus, we investigate diagrams which have only one <^u in the external lines. 

At first, we focus on marginal operators in the effective potential. The most dangerous 
marginal operator for the instability of the VEV hierarchy is Agl^^K^^^;^) (+h.c.), which is 
induced from diagrams in Fig.l (a) and (b). It is because this operator breaks Z2-parity and 
induces linear term of f 2, which might possibly destroy the VEV hierarchy. Here we note that 
Fig.l (a) and (b) are only 1-loop diagrams which induce Agl^^K^^^;/) (+h.c.). Neither lepton 
nor quark 1-loop diagrams contribute Ag due to the Z2-parity, since one additional external ^^ 
needs one additional right-handed neutrino propagator inside a loop which requires one more 
(^y. Fig. 1 (c) and (d) induce another Z2-parity violating operator, A7|$^|($'''$y) (+h.c.). 
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Fig.l: Z2-violating 1-loop diagrams. 
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Figure 1 (a) ((c)) is expected to dominate (b) ((d)) because of |Ajp :» (725 so Ae and A7 are 
estimated as 



3A1A5 



m% 



47]-2 
3A2A5 



m^ 



ni 



2\2 



?7^^ 



m\ + ?7i2 In 



mt 



\m% 



ni% 



47r^ (m2 — fn\ 



2\2 



2 2 I 2 1 

1712 "~ '^l + '^l ^^ 



mt 



m? 



(2.28) 
(2.29) 



Taking into account all irrelevant operators which have only one (f)y in the effective operator, 



correction for \W might be of order 



3 

27r2 



log 



HI 

V2 I 



at most. This correction contributes the 



stationary condition of V2 in Eq.(2.4), and modifies it as 



6„,3 



= m2f2 — ra^vi + A2t'2 + ^v^V2 H vi + 



3A 



7 2 



(2.30) 



Remind again that tiny VEV of |f2|(^ K'i|) is originated from tiny term of m\vi. Thus, an 
induced term from the radiative correction of ^v\ must be smaller than m\vi to preserve the 
VEV hierarchy. Actually, the ratio of them is estimated as 



XqvI 



2m^vi 



3 

4^ 



log 



V2 



(2.31) 



at most. This means that the order of |f2| is not changed but its factor might be modified about 
0.8 (2) by the radiative corrections in Majorana (Dirac) neutrino scenario. This magnitude 
comes from a maximal (may be over-) estimation, and anyhow, the orders of VEVs are not 
changed. (Actually, this modification becomes much smaller about (9(1)%, if we use Higgs 
self-couplings of (9(0.1).) Thus, the VEV hierarchy itself is stable against radiative corrections. 
As for higher-loop effects, they are at least suppressed by an additional loop-factor j^, and 
we cannot find any diagrams which have larger contribution than above diagrams. Therefore, 
the VEV hierarchy itself is stable against radiative corrections, and we can conclude radiative 
corrections do not destroy the VEV hierarchy in both Dirac and Majorana neutrino scenarios. 



3 SUSY neutrinophilic Higgs doublet model 

In this section, we analyze vacuum structures of Higgs potential in the SUSY neutrinophilic 
Higgs doublet model at first, and next investigate a stability of this VEV hierarchy against 
radiative corrections. 



3.1 Vacuum structure in tree-level potential 

The SUSY neutrinophilic Higgs doublet model has four Higgs doublets[l2l[T3], and the super- 
potential is given by 



W 



y'^Q'-HJjR + y^QLHdDn + LHdER + y^LH.N 



(3.32) 



where H^, gives Dirac neutrino masses and H,^/ does not couple with any matters. Note 
that Hu and Hd are usual MSSM Higgs doublets. This superpotential is for Dirac neutrino 
scenario, and Majorana neutrino scenario can be realized when Majorana mass of right-handed 



neutrinos MN"^ is included in Eq.(3.32). The Z2-parity assignment of the fields is shown in 
the following table. 



Fields 


Z2-parity 


Lepton number 


MSSM Higgs doublets H^^Ha 


+ 





neutrinophilic Higgs doublets H^, Hyi 


— 





right-handed neutrino A^ 


— 


1 


others 


+ 


± Ideptons, 0:quarks 



Note that Z2-parity is softly broken by p, p' , where |p|, \p'\ <^ \p\, \p'\. The Higgs potential is 
given by 

V ^ {\p\' + \p?)HlH^ + {\p\' + \p'\')H\Hd + (|/.f + \p'?)HlH, + (l/xf + \p?)hIh^, 

+ f {Hi\H^ - H\\Hd + HI\h^ - Hl^H^^j 

, ^2 / _a ^a ^a _a 

a ^ 

—m 



+BpHu ■ Hd + B'p'H, ■ H,, + BpH^ ■ H,, + B'p'H, ■ H^ 
+p*pH\H,, + p*p'HiH, + p'*p'Hl,Hd + p*pHlHu + h.c. 



(3.33) 



where r" and dot mean a generator and cross product of SU(2), respectively. B'p' {Bp) cor- 
responds to m\ in the non-SUSY z/THDM, and its smallness plays a crucial role of generating 



tiny VEVs oiH^y. The magnitude of \B'p'\ {\Bp\) is (C(10~°-5) GeV)^ for Majorana neutrino 
scenario, and is ((9(10^^) MeV)^ for Dirac neutrino scenario. We assume {—m'jj^) < for the 
suitable electroweak symmetry breaking and real VEVs as 



(H.,) 



«.w^';;.w^ «!.(//.) 







(3.34) 



By taking /i, p, 5-parameters to be real and denoting M^ = |/ip + |pp — m'jj^{< 0), Mj 



/|2 



m 



stationary conditions are given by 



Hd 



;> 0), M,2 ^ \^f + \pf-ml(> 0), and M^, ^ |^f + |p|2 + ^2 (> g)^ ^^e 



1 f~i\ / 1 

0=2^ = ^«^" + 4(^? + ^2)^«(^« -'"d + '^l- '"I') - ^/^^rf - Bpv^> + ifip' + p'p)v^, 
0=2^ = M^Vd - ^{9i+ 9l)vd{vl -vl + vl- vl,) - Bpvu - B'p'v^ + {pp + p'p')v^>, 
= -^ = MX + 4(^1 + ^2)^.(^« - ^d + ^' - ^'') - ^W - ^'P'^d + if^p' + p'p)vu, 
= 2^^ = ^^'^^' - 4(^1 + ^2)^!-'(^« -'^d + '^l- ^l') - B'p'v^ - Bpvu + {pp + p'p')^rf- 



Let us investigate the vacuum structure with a parametrization of Vu = f sin/3cos7,t'd = 
vcos/3cos7,t;jy = wsin/3i,sin7,f,y/ = f cos /S,^ sin 7. At first, we focus on the vacuum which 
neutrinophilic Higgs doublet model requires, i.e., \vu\, \vd\ ^ \vu\, \vu'\- This condition induces 
the usual MSSM relations for Vu.Vd as 

Ml - -^{gl + giy cos 2/3 - Bp cot /3 ~ 0, M^ + -{gl + gW cos 2/3 - Bp tan /3 ~ 0, 



which means 



^2 ^ 



^? + ^i V cos 2/3 
They induce tiny v^,v^' through tiny p, p' as 



A/f2 — M^ 

" ^-(M„^ + M,^)), sin 2/3 



2Bp 



M^ + M]- 



(3.35) 



[M2 - lisl + gl){vl - t;2)] B'p'vd - {pp' + p'p)^^^ 



+ B'p' 



BpVu - {pp + p'p')vd 



[Mi + \{gl + gl){vl - vl)] [Ml, - \{gl + gl){vl - vf)] - B'^p'^ 



[Ml + \{gl + gl){vl - vD] Bpv^ - {pp + p'p')vd 



v,/ 



+ B'p' 



(3.36) 
B'p'vd - {pp' + p'p)vu 



[Ml + \{gl + gl){vl - vl)] [Ml, - \{gl + gl){vl - vD] - B'^p'^ 



(3.37) 



At this vacuum, the potential height is estimated as 



V ~ v\Mt, sin^ /3 + M| cos^ (3 - 25/i cos (3 sin (3) + -{gi + gi)v* cos^ 2/3. 



(3.38) 



Next, we study the vacuum at |w„|, \vd\ ~ \vu\, \vi,i\. Where, by neglecting both p and p', 
the stationary conditions become 



^l-Vu + ^(^1 + 9l)vu{vl -vl + vl-vl,)- Bpvd = 0, 

MdVd - -^{gl + gl)vd{vl -vl + vl- vl,) - BfiVu = 0, 

MX + \i9! + gl)vu{vl -vl + vl- vl) - B'p'v,, = 0, 

Ml>Vu' - ^{gl + gl)vu'{vl -vl + vl-vl,)- B'p'vy = 0. 



(3.39) 
(3.40) 
(3.41) 
(3.42) 



It is easy to show that only v^ = v^ = can satisfy the stationary conditions in D-flat direction 
of v^ = vl- 

Numerical analyzes show that the vacuum aX v,^ = v^ = is the global minimum in suitable 
parameter regions of \B'\, |p'| = 0(10^) GeV and positive M^,M^> = O{10^) GeV. This resuh 
is originated from an initial setup that only soft mass squared of H^ is negative. (See, case (1) 
of Table 1 in z/THDM.) Similarly, we can show that there is no vacuum at \vu\, \vd\ ^ \vu\, \vu'\- 
Anyhow, the vacuum exits only at |f„|, \vd\ ^ |fjy|, \vu'\, which is the desirable vacuum in the 
neutrinophilic Higgs doublet model. 

3.2 Stability against radiative corrections 

Let us investigate the stability of the VEV hierarchy against radiative corrections in the SUSY 
neutrinophilic Higgs doublet model. As in non-SUSY case, we can estimate 1-loop radiative 
corrections in a SUSY effective potential. 
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Fig.2: Z2-violating 1-loop diagrams in SUSY. 
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The most dangerous marginal operator in the effective potential is \'{H'lHy){H\Hd) (+h.c.), 
which is induced from Z2-violating diagrams in Figs. 2 (a)~(c). The absolute value of A' is 
roughly estimated as ^|^^| at most, where ?7i is a Higgs mass in a loop. Notice that neither 
(s)lepton nor (s) quark contribute A' at 1-loop level due to the Z2-parity similarly in non-SUSY 
i^THDM. It is because one additional external H^ needs one additional right-handed neutrino 
propagator inside a loop, which requires one more external H^. Anyhow, this term modifies 



the stationary condition of Vy in Eq.(3.35) as 



Q = Mlv,--{gl + gl)v, 



2 2,2 2 I 



2A' 



VuVl 



{91 + 92) ^^ 



- B'p'vd + ipp' + p'pM3A3) 



Taking into account all irrelevant operators which have only one H^, in the effective operator. 



correction for |A'| might be of order £^ 



92 \B'p' 



log 



'"u,d I 



at most. Remind that tiny VEV of v^, is 



originated from the small mass parameters of B' p' as in Eq.(3.36). Thus, in order to preserve 



the VEV hierarchy, ly^^wf^l must be smaller than \B'p'vd\ in Eq.(3.43). And, this ratio is 
estimated as 



(3.44) 



\/ 2 

2B'p'vd 


9i 

87r2 


VuVd 


log 





This value is too small to influence the stationary conditions in both Dirac and Majorana 
neutrino scenarios. We can also show that higher-loop diagrams induce smaller corrections 
due to the loop suppression factors. Therefore, we can conclude that the potential is stable 
against radiative corrections in SUSY neutrinophilic Higgs doublet model. 

4 Summary 

A neutrinophilic Higgs model has tiny VEV, which can naturally explain tiny masses of neu- 
trinos. There is a large energy scale hierarchy between a VEV of the neutrinophilic Higgs 
doublet and that of usual SM-like Higgs doublet. In this paper, we have analyzed vacuum 
structures of Higgs potential in both SUSY and non-SUSY neutrinophilic Higgs models, and 
next investigated a stability of this VEV hierarchy against radiative corrections. We have 
shown that the VEV hierarchy is stable against radiative corrections in both Dirac neutrino 
and Majorana neutrino scenarios in both SUSY and non-SUSY neutrinophilic Higgs doublet 
models. 



Note added 

After preparing our submission of this paper, we notice a paper [19], where authors also 
analyzed the vacuum stability against radiative corrections in the non-SUSY z/THDM with 
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Dirac neutrino scenario. Their results are consistent with ours. They calculated the 1-loop 
effective potential and the quantum corrections to VEV hierarchy. On the other hand, we 
estimated the most dangerous contributions to the VEV hierarchy and confirmed the stability 
also in SUSY and Majorana cases. 
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